Quantum dynamics in a time-dependent cylindrical trap 
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Solutions of the Schrodinger equation is examined for a particle inside a cylindrical trap of a 
circular time-dependent cross-section: the analytical expressions for the energy and momentum 
expectation values are derived with respect to the exact solutions, and the adiabatic and the sudden 
change of the boundary are discussed. The density profile in the time, in a given observation 
point, resembles the diffraction-in-time pattern observed in a suddenly released particle but, with an 
enhanced fringe visibility. Numerical computations are presented for both contracting and expanding 
boxes. 
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^vq , Quantum mechanical particle motions which are subject to time-dependent boundary conditions belong to a main 
class of the time-dependent dynamical systems. Finding the exact solutions of the Schrodinger equation in these 
time-dependent systems, is not an easy task [Tj-Q- 

Diffraction in time [lj], quantum temporal oscillations of matter waves released from a confining region, is just one 
of the very interesting effects that are seen in such systems . For a beam of particles which is incident from the left 
on a moving wall the visibility of the fringes is enhanced in comparison with the motionless case @ . 

Moshinsky's theoretical work has been studied for particles which are suddenly released from a ID box from a 
^ | spherical Q and from a cylindrical trap @. 

Exact solutions of the Schrodinger equation for a particle in a ID box with a moving wall has been found @, H|. 
Using the semiclassical approximation, Luz and Cheng (lp| evaluated the exact propagator of the problem. 

It has been noted by Godoy and Okamura Q, that particles motion in a nanoscopic circuit can be regarded as a 
quantum wave moving inside a wave guide filled with a scattering medium, and transient currents in these circuits can 
be produced by rapidly removing (or adding) boundary walls to the guide. See [Tl| for a recent review on quantum 
transients. 

With this in mind, we want to examine the solutions of the Schrodinger equation for a particle inside a circle with 
a varying radius. There are two extra motivations for this study. First, the sudden removal of the boundary is an 
C"-"- idealized case. Second, possible applications to optical effects connected with the moving mirrors are additionally 
encouraging. Although, the limit of infinite velocity of the moving boundary is clearly corresponds to the sudden 
removal of the wall, due to the non-relativistic nature of the Schrodinger equation, the speed of the wall has to be 
very less than the speed of light, u <C c. 
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II. EXACT SOLUTION 

Consider a particle with mass p inside a cylindrical wane guide of circular cross-section. The cylinder has a radius 
L{t) = a + lit, which uniformly changes with time, and its axis is taken along the z-dircction. The potential energy 
function is zero if r < L(t) and infinite otherwise. The 2D Schrodinger equation for the transverse wavefunction in 
the polar coordinates r = (p, <f>) is then 

, , d . . h 2 n d ( d \ Id 2 



p dp \^ dp ) ' p 2 dej) 2 



*(r,t) . (1) 
In the absence of the boundary motion, u = 0, the solutions of eq.([T]) would be 
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with eigenvalues E mn — h 2 x 2 nn /(2pa 2 ), where x mn are the roots of Bcssel functions, Jm{x mn ) — with m = 0, 1, 2, ... 
and n = 1,2,3,.... All Bessel functions with ra^O have a zero at the origin, but one has not to consider these zeros 
to have a non-zero wavefunction. Thus, the instantaneous energy eigcnfunctions and eigenvalues would, respectively, 
be 



l 



L(t) \Jm+l(Xmn)\ 



im(j> 



Lit)) 



2fx 



L 2 (t) 



(3) 
(4) 



which satisfy the relation 



ih Ft { 



if* 

-i-^ J dt E mn (t) 



u mn (r,t) = E mn {t) cxp 



-i- I dt' E mn (t') 



u mn . k (r,t) 



Exact solutions of the Schrodinger equation ([T]) for the above problem have been found [l2j , 



*rrmO,<M) = C 



L(t) 



exp 



iH p 
■u 



2Ti L(t) 2p rnn aL(t) 



(5) 



(6) 



Unknown coefficient C is determined by the normalization condition J Q dpp J d<p\^ mn (r, t)\ 2 = 1, which apart 
from a constant phase factor yields in C — s/2/J m+ i(x mn ), by using the orthogonality of the Bessel functions fl3| 

ds S J m (x m pS)t/ m (x m gS) — [t/j^-j-l (^mp)] &pq ■ (7) 



Thus, exact solutions read 

* mn (p,0,*) = exp 



L(t) 



4a 



(8) 



where we have introduced dimensionless quantities a = p,au/{2%) and £(£) = L(t)/a. It must be mentioned that 
including the motion in the z-direction, would introduce an extra factor exp[ifcz — ihk 2 t/(2p)}. Functions \E' mra (r, t) 
vanish at p = L(t), remain normalized as the boundary moves and form a complete orthogonal set. Thus, the general 
solution of eq. ([IJ can be expanded in terms of functions (f8|). 
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with time-independent coefficients c m i n i determined from the relation 
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General solution can also be written as a superposition of instantaneous eigcnfunctions, 
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*(r,t) = rfm '™' (*) exp / dt' E m , n ,(t') u m > n >(T,t)= ^2 ^Z b m'n'{t)u m ' n '{r,t) 

m'—0 n' — l m'-On'-l 
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now with time-dependent coefficients b m f n * (t) determined from the relation 
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L(t) 



2tt 



dpp d<t>U* m , n ,{Y,t)^{Y,t) 



Using eqs. ^ and (|12p and the orthonormality of functions exp (imcf>) /y/2n, one obtains 
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where 



f 1 2 

>(t,Ct)= ds S e - M «(') s Jtn'(x m ' n 's)J m '(x m ' n "s) , (14) 



with the property I mnn >{t, -a) = I mnn ,(t, a). 

If the particle is initially in an energy eigenstatc, i.e., ^(r, 0) = u mn (r, 0), then 



8mm' TJ 7 v~ j / T7^mnn'(0, a) , (15) 



which is not an unexpected result as the angular momentum operator L z — —ih d/d(f> commutes with the Hamiltonian 
of the system. For a«L one has l mnn >(0,a) ~ 8 nn i(J m+ i(x mn )) 2 /2, and thus from eq. (|9|), one obtains ^(r, i) ~ 
f mn (r,f), which is in agreement with the result of the adiabatic approximation. 

For completeness we should add the relations of the energy expectation value and the propagator of the problem. 
The expectation value of the energy of the particle is obtained from 

,L(t) ,2tt q 

(#>(*) = J dppj^ d<f>**(r,t)ih-*(r,t) 

= J2 Mm'n'WI 2 E m , n .{t) = J2 \ b ™>n>(t)\ 2 E m , n ,{t) . (16) 
rn'n' rn'71' 

Propagator of the problem can be constructed as follows: 

i*(t)> = s(t,toi*(to> = EEi*-»( i ))^-™wi 5 (*>* / )i*-'»'(t'))(*m' n '(t / )i*(* / )) 

ran m'n' 

mn 

where S(t,t') is the time evolution operator and we have used the fact that if the particle is in the sate \^ m n) a t t' , 
it remains in that state as the wall moves, i.e., S(t,t')\ty mn (t')) = \^ m n(t)). Now, we write this equation in the form 

*(r, t) = f 1 dp'p' r dcb'K(r, t- r', t')*(r', f ) , (17) 
Jo Ja 

where we have introduced the propagator as 

oo oo 

K (r, t; r', t') = E E *™"( r ' *)*™»( r '> 
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L(t)L(tf) ^ [J^+lC^n)] 2 

I/km / p 2 p' 2 \ ih x mr ft t' 



2k \L(t) L{t') 2p a \L(t) L{?) 



x exp 

Expectation value of an arbitrary operator A is given by 

oo oo oo oo 

(*|i|*>(t) = EEEE<.vU*™v(#4»- (19) 

m'—Q n' — l m—0 n— 1 

where 

(* m '„/(t)|A|* m „(t)) = / #/ dp P W m , n ,(p,<l>,t) A(p,</>) ^ mn {p,4>,t) , (20) 

Jo Jo 

are the matrix elements of the operator with respect to the states ^> mn (r,t). It must be noted that the expectation 
values can also be found by expanding the wavefunction in terms of the instantaneous energy eigenfunctions as we did 
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in eq. (|16p . Our aim is the examination of the uncertainty relation. To this end, we at first compute matrix elements 
of the position, the momentum and the Hamiltonian operators. 

The position and the momentum operators in the arbitrary radial direction po = x cos (f>Q + y sin 4>q , arc 



q = xcoscpo+y sin (p , 
P a = Px cos 0o + p y sin 4> Q 



COS (pQ COS 



dp 



1 • A Q 

-sm^— 

p dq> 



sm 0o sin 



dp 



1 ^ 5 



which are canonical conjugate variables, [q ,p ] = ih. By straightforward algebra, one obtains 



(*tn'n'(t)l9ol*mn(t)> = 
(*m'n'(*)|j»ol*mn(*)> = 



(*™'T.'(t)l9*l*mn(t)» 



\ ^ r 



,(t)\p 2 \* mn (t)) 
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4a 2 A (3) , 

-"-mn'n 



1 



■(m 2 A 



2 4 (-D 

mn'n 



D 



(0) 



£ 2 W 

rf/ 5 P *m'n'(/°i^*) 



2^ Jo 

1 (tt m ' n '(t)|pg|tt mn (t)) 
7T 2/i 



1 5 



mn'n' 

90 2 _ 



^ (1 



) 

mn'n 



D 



(2) 



p dp \ dp 



8 



(21) 



(22) 
(23) 

(24) 



(25) 



(26) 
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(fc) 



C! 



(fc) 



(fc) _ 



ds S J m (yXjYifi' s) Jmi^mnS) 



ds S Jmi^mn' s) , Jm^mn$) ) 

as 

/ (is S Jmi^mn' s) Jmi^mn^) ; 

Jo «« 



(27) 
(28) 
(29) 



are constant coefficients. As one sees, all of the expectation values are independent of </>o, the direction we started 
with it. A^ n , n is symmetric with respect to two indices n! and n, while B^ n , and C^ n , n do not have any symmetry. 
It must be mentioned that A^ m ^) n diverges for m = 0, but it is not problematic because the wave- function is not 
-dependent for m — and thus, the term m 2 A m }) n will not appear in such a case. 



From cqs. (|19[) . (|22[) and ([23]), one sees that position and momentum expectation values are zero, irrespective of 
the shape of the wavefunction, (^\q \^)(t) = and (*|p |*)(<) = 0; and from eqs. ([19} and ([26|). (*|iJ|*)(t) = 
(y\p 2 o W(t)/(2np). 



III. UNCERTAINTY RELATIONS 



Let us suppose that one succeed to construct the initial wavefunction to be the state \D' mTI (r, 0). Such a task 
can be done by appropriately superposing instantaneous energy eigenstatcs. If the particle is initially in the state 
*(r,0) = *mn(r,0), then 



C-m'n' 
^m 'n' (J*) 



$ mm' $nn' j 

8* 



Jm-j-l {Xmn) J m+1 i^mn' ) | 



exp 



1:1 - 



4a 



(30) 
(31) 



From eqs. © and (|3T)|) . one finds that if ^(r,0) = 4' m „(r, 0), then time-evolved wavefunction would be $ mn (r,t) in 
any instant of time. 
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By using eqs. ([TBI and (j3"Tj) . one obtains 

(H) mn (t) 



L(t) 



'(0,a)| 5 



(32) 



for the energy expectation value. From this relation it is not obvious how the energy expectation value changes with 
time, whereas eq. (|34[) is completely informative in this connection. 
Using eqs. (|22l), ([23]l. ([24)) and (f25]). one obtains 



(q 2 ) mn (t) 
(p?L»(*) 



2irA 



(3) 



2- 



2 ' 



a 2 |J m +i(a; m „)| 2 



4a 2 A (3) + — !— fm 2 ^" 1 ) -C (1) ) 

^ mnn 1 <~2 ^ mnn ^mnn > 



em 



(2irn)(H) mn (t) , 



A<? , ra „(i) = ^<g 2 L„W -<<z„>L(*) = 



2tt4 (3) 



m+1 



(^mn ) | 



A PoimB (t) = V^L»(*)-<Po>L(*) = 

A9 , mn (t)Ap„ ran (t) 



Ua 2 



Ag , m „(t) 
L(t) 



1 2 7 r(m 2 AL! ) 



G, 



(i) 



£ 2 (£) |J m +l (£„„)! 2 



47T 



2 | t/m-fl (^mn ) | 2 



(33) 
(34) 

(35) 

(36) 
(37) 



0. 



where, for abbreviation expectation values have been shown as (...) mn (t) instead of (\P mn (f)|...|\D r mn (t)); and Bmnn 
In the appendix the solutions of Amnn and Cmnn will be given in terms of the hypergeometric and the regularized 
hypcrgcometric functions. 

Considering the above equation, one finds 

i) Aq mn (t) changes linearly with time, it increases for expansion while decreases for contraction. This is an expected 
result. 

ii) There is an additional time- and ^-independent contribution to Ap a mn in comparison with the stationary bound- 
ary. The second term under the radical sign is time-dependent and decreases with expansion while increases in 
contraction. For an expanding circle, Ap mn (i) becomes constant in the limit t — > oo. 

iii) Eq. (|37[) shows that the product of the uncertainties is time-dependent, a result which is absent in 
stationary systems. According to the general uncertainty principle Aq,Api > h/2, where qi and Pi 
arc the non-commuting i— component of the conjugate variables r and p. This means the quantity 



47tv/ (m 2 A 



(-i) 



/| J m +i{x mn )\ 2 , which is the product of uncertainties in the corresponding stationary 



I (3) 



system, must be greater than 1. 



IV. DISCUSSION 



As eq. ([3~4")) shows and fig. [JJ confirms, expectation value of the energy decreases with time for an expanding box 
while increases for a contracting one, as time elapses. This is an understandable result in the context of the old 
quantum theory [l^ , as is reasonable by considering the uncertainty relations ([33]) and (|3"6"|) fl(| . 

Fig. [2] shows that, in the presented region, Amnn/\Jm+i(x mn )\ 2 and A^ nn /\J m+ i{x mn )\ 2 becomes constant for 
n > 10, while |Cmnn|/| Jm+i{xmn)\ 2 increases with n. In addition, A^mnn / \Jm+i{x mn )\ 2 is approximately independent 
of quantum number m for n > 10, while Cmnn/\J m +i(x mn )\ 2 is approximately the same for all values of n and 
dominates for n > 3. From these considerations; and eqs. (|33p and (f3~4"|). one finds that at a given time and a given 
wall velocity, (q 2 ) mn becomes independent of n and m for n > 10, while {p 2 ) mn is independent of m, but increases 
with n. 

For numerical calculations of the wavefunction, we define two new quantities X mn = 2na/x mn and v mn = E mn / (2-kTi) 
and from them new dimensionlcss position and time coordinates r\ — p/X m n and T = v mrl t\ and dimensionlcss radial 
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FIG. 1: (Color online) {H) ol (t)/{H) ol (0) versus £(t), for (a) contracting and b) expanding circular boxes. Expectation value 
of the energy of the confined particle increases (decreases) for a contracting (an expanding) box. 



probability density g m n(ViT) = A^^i?^, T)\ 2 , where R(j],T) stands for the radial part of the wavefunction. With 
these quantities, the location of the wall is determined from the relation 

a 

£(T) = 1 + 27T— T , 
or 

where, 

[MX fid flXnm %mn 

2n 2n ixa 2 

corresponds to the wall velocity that is equal to the initial classical velocity of the confined particle. 

We have plotted g mn (r],T) in fig. [3] for a particle initially in the state wo,io, against r\ at different times. 

A system undergoes an adiabatic evolution when t\ <t c , where t\ (t e ) is the time-scale over which internal (external) 
variables of the system changes [13] • For our particle-in-a-box, t e = a/u and = a/v mn . Figures |4] and [5] show that 
the particle, initially in an instantaneous energy eigenstate will remain in the same eigenstate when |a| <C a mn , 
as a consequence of the adiabatic approximation. Figure [4] shows that sudden approximation does not work for a 
contracting well [l5[ • 

Figure [6] represents density probability versus time in a given observation point. 

T\ and T-i are dimensionless classical flight times from the front and back edges of the circle to the dimensionless 
observation point. In contrast to classical mechanics, when a > a mn , one sees a non- monotonous increasing behavior 
of the density for T < 7\. A non- monotonous decreasing behavior is seen for T > T 2 , irrespective of the wall 
velocity. Difference between the adjacent cxtrcmums (visibility) decreases with a. The constructive interference with 
the reflected components from the moving boundary for a < a mn , increases the fringe visibility. At long times, the 
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FIG. 2: (Color online) (a) A\nnl/\J m +i{x mn )\ 2 , (b) A\Z' nn /\J m+1 (x mn )\ 2 and (c) |C^A„|/| J m +i(x mn )\ 2 versus quantumnumber 



n for different values of quantum number m. 



behaviour of the density in the observation point, is the same for all values of the velocity parameter a. This can 
be understood by considering the behavior of functions ^ m n at long times: at this limit at an arbitrary observation 
point r = r , the coordinate-deepened phase goes to zero. Thus, from eqs. © and (JTUJ), one obtains 



lim \ip(r ,t)\ 



which ultimately becomes zero. 



^^o,^) 



ut I J, 



m+l 



(x m n)V Jm{Xmn ut ] 



(38) 



SUMMARY 



Exact solutions of the Schrddingcr equation for a particle in a circular hard box with a moving wall in uniform 



motion, contain a coordinate-dependent phase exp 



' 2h L(t) 



which is absent in the corresponding stationary 



boundary problem, whose significance has been already emphasised [1, [T3| in ID systems. Propagator of the problem 
is constructed and the matrix elements of the position, the momentum and the energy observables were derived 
with respect to the exact solutions. It was seen that the uncertainty product increases with time due to a timc- 
and velocity-dependent term; and expectation value of the energy increases (decreases) with time, for a contracting 
(an expanding) boxes which is consistent with the uncertainty relations. The density profile in the time, in a given 
location, resembles the diffraction-in-time pattern observed in a suddenly released particle, but it actually shows an 
enhancement of the visibility of the fringes, the 2D version of the EDIT (Enhanced diffraction in time) reported in 
ID in We have provided a full characterization of the quantum transients in an expanding/contracting cylindrical 
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Dimensionless position 



FIG. 3: (Color online) Dimensionless radial probability density 00,10 (*7, T) for a — ao,io versus dimensionless position coordinate 
r\ at the instant of time when the wall arrives at (a) £ = 1, (b) £ = 1.5, (c) £ = 2, (d) £ = 2.5, (e) £ = 3 and (f) £ = 3.5. 



box resulting from the breakdown of adiabaticity. We close by noticing that, should it be desirable, it is possible to 
suppress quantum transients in this system by engineering a shortcut to adiabaticity fl8j . 
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VI. APPENDIX 

Here, we write integrals Amnl , AmL and Cmnn in terms of generalized and regularized generalized hypergeometric 
functions: 

1 f 1 

2m Jo 

= 4- m (2m-l)! 2 F 3 (m,m+^m+l,rn+l,2m+l-,-(x mn )^j if m ^ , (39) 

(3) f4- m m(m + l)(a;„ m ) 2 (2m- 1)! 2 f 3 (m + i, m + 2; m + 1, m + 3, 2ra + 1; - (s mn ) 2 ) if m + 0, 

1 i 2 ^ 3 (|,2; 1, 1, 3; -(x „) 2 ) if m = 0, 1 Uj 
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FIG. 4: (Color online) Dimensionless radial probability density Qo,i(ViT) against dimensionless position coordinate rj at the 
instant of time when the wall arrives at £ = 0.1, for six different values of contraction rate; a — —0.01 ao,i (black curve), 
a — — ao.i (red curve), a = —5 ao,i (green curve), a = —10 ao,i (blue curve), a = —15 ao,i (yellow curve) and a — —20 Qo,i 
(magenta curve). For a rapidly moving boundary, the sudden approximation does not work for a contracting box, but, in the 
opposite limit the process is adiabatic. 



^ hi n i 



and, 



ds s 



ds 



(jEmn ) 



2 r 



(x mn ) 2m 3F4 (m + g,m + l,m + l;m,m + 2,m + 2, 2m + 1; -(x mn ) 2 ) 

4 m (m - 1)! m! 

24-2(^)4-1(1™™) (4m(m 2 - 1) + (x„ m ) 2 (l - 2m)) 2 (x m „) 2 - 2m(m + 1) 



+ (Jm-2(x m „))' < 



"4" («/m — 1 (^mn)) I 



1 4m(m- l)(2m 2 - 2 - (a; m „) 2 )\ (J m+ i(x m „)) 



(*^mn ) 



(^mn) 
if 777, 7^ 



(41) 



(7(1) 



:a;on(J'i(a;on)) 2 if m = , 



(42) 



where the generalized hypergeometric function p Fq(a\, a p \ b\, 6 9 ; z) in terms of the Pochhummer symbol |13j . 

(a + n - 1)1 



(a), 



(a -1)1 



, (a)o = 1, 
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Dimensionless position 



FIG. 5: (Color online) Dimensionless radial probability density £0,2(77, T) against dimensionless position coordinate rj at the 
instant of time when the wall arrives at £ = 2, for five different values of expansion rate; a = (black curve), a = 0.01 Qo,2 
(red curve), a = ao,2 (green curve), ct = 5 ao,2 (blue curve) and a = 10 Qo,2 (magenta curve). Both the adiabatic and the 
sudden approximation work for an expanding box. 



becomes 

(ai)fc ■■■ (a p )k z h 



p F q {ai,...,a p ;bi,...,b q ;z) = ^ 



^ (h) k ...(ft,)* fc! ' 

and the regularized generalized hypergeometric function p F q (ai, a p ; &i, b q ; z) is 

pF q (ai, a p ; 6i, b q ; z) 



p F q (ai, ...,a p ;b lt ...,b q ; z) 



r(6i)...r(6,) 



and r(z) = J °° t z 1 e '(if is the Gamma function. Recurrence relations of the Bessel functions have been used in (|39l 
to re-write the integrand, which is more appropriate for numerical calculations. 
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Dimensionless time 



FIG. 6: (Color online) Dimensionless radial probability density Q mn (ri,T) for a particle initially in the state (a) Mo, 6 and (b) 
Mo, 15) against dimensionless time coordinate T at dimensionless observation point n/n, for three different values of 

velocity parameter; a — 0.9 a mn (black curve), ct = a mn (red curve), a = 2 a mn (green curve). T\ and T2 are dimensionless 
classical flight times from the front and back edges of the circle to the dimensionless observation point. As the velocity of the 
moving boundary increases, the fringe visibility reduces. 
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